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Q_i| We perform a detailed study of the process e + e — > ^V[Vi and its sensitivity 

to anomalous gauge boson couplings of the ^yWW vertex. We concentrate on LEP 



II energies, \fs = 200 GeV, and energies appropriate to the proposed Next Linear 
Collider (NLC) high energy e + e~ collider with center of mass energies ^fs = 500 
and 1 TeV. At 200 GeV, the process offers, at best, a consistency check of other 
processes being considered at LEP200. At 500 GeV, the parameters k 7 and A 7 can 
be measured to about ±0.05 and ±0.1 respectively at 95% C.L. while at 1 TeV, they 
can be measured to about ±0.02. At the high luminosities anticipated at high energy 
linear colliders precision measurements are likely to be limited by systematic rather 
than statistical errors. 
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I. INTRODUCTION 



The major preoccupation of particle physics is the search for physics beyond the stan- 
dard model or equivalently, for deviations from standard model predictions. To this end, 
measurements at the CERN LEP-100 e + e~ collider and the SLAC SLC e + e~ collider |l| 
have provided stringent tests |||| of the standard model of the electroweak interactions 
Q. However, it is mainly the fermion-gauge boson couplings that have been tested and 
the gauge sector of the standard model remains poorly constrained. A stringent test of the 
gauge structure of the standard model is provided by the tri-linear gauge vertices (TGV's); 
the ■yWW and ZWW vertices. Within the standard model, these couplings are uniquely 
determined by SU(2)l x U(l) gauge symmetry so that a precise measurement of the vertex 
poses a severe test of the gauge structure of the theory. If these couplings were observed 
to have different values than their standard model values, it would indicate the need for 
physics beyond the standard model. 

The study of the trilinear gauge boson couplings by studying W pair production is one 
of the primary motivations for the LEP200 upgrade [0-0] with a precision of 30-40% is 
expected from cross section and W angular distribution measurements. In the far future 
there is growing interest in the physics that can be done at high energy e + e~ colliders with 
\fs = 500 GeV or s/s = 1 TeV, referred to as the Next Linear Collider (NLC), the Japan 
Linear Collider (JLC) or the CERN Linear Collider (CLIC) @-|13[- Various options are 
being studied including and 77 collisions where the energetic photons are obtained by 
backscattering a laser on one of the incident leptons. Measurements at these colliders are 
very sensitive to anomalous couplings with e'j and 77 collisions putting some of the more 
stringent bounds on anomalous WW'j couplings |14|-[16|| . 

A problem common to many processes used to study TGV's is that they involve both 
the WW'y and WW Z vertices making it difficult to disentangle the contributions. In a 
previous paper we presented a detailed study of the process e + e~ — > viVi\i r fi~ motivated 
by our interest in isolating the WW Z and WWj vertices by appropriate kinematic cuts 
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on the invariant mass of the fi + [i~ W% - Included in this final state are contributions from 
the underlying process e + e~ — > viv^* — > b>ii?ifi + fi~ which shows up most dramatically when 
Mfi+fj,- — > 0. However, because of the muons' masses it does not quite isolate the process 
we are interested in. In this paper we take the obvious limit and study the sensitivity of the 



process e + e — > to anomalous WW'j couplings H,|TH]- This process has also been used 



as a means of counting the number of light neutrino species [19 



To parametrize the WW^ vertex we use the most general parametrization possible that 



respects Lorentz invariance, electromagnetic gauge invariance and CP invariance P, |20|J21 
since it has become the standard parametrization used in phenomenology and therefore 
makes the comparison of the sensitivity of different measurements to the TGV's straightfor- 
ward. We do not consider CP violating operators in this paper as they are tightly constrained 
by measurement of the neutron electron dipole moment which constrains the two CP vi- 
olating parameters to |/c 7 |, |A 7 | < (9(10~ 4 ) |22fl . Therefore the WlVj vertex has two free 
independent parameters, k 7 and A 7 and is given by [fj pOfl : 

Cww, = -ie {(W^WM* - W^W^) + k^W^ - A-W^W^j (1) 

where and represents the photon and W~ fields respectively, W^ v = £? M W„— dJW^ and 
Fp,v — d^A u — d u A^ where A is the photon and M\y is the W boson mass. Higher dimension 
operators would correspond to momentum dependence in the form factors which we ignore. 
At tree level the standard model requires k 7 = 1 and A 7 = 0. Note that the presence of the 
W-boson mass factor in the Ay term is ad hoc and one could argue that the scale A of new 
physics would be more appropriate. We will conform to the usual parametrization and will 
not address this issue any further. 

We studied the sensitivity of this process at TRISTAN and LEP/SLC energies where 



there exists data |23| , |24|| that we could in principle use to bound the WW^j couplings. How- 
ever, we found that the process was insufficiently sensitive at these energies to put meaningful 
bounds on the WW'j coupling with the integrated luminosities already accumulated or ex- 
pected in the foreseeable future. We therefore start with y/s = 200 GeV appropriate to 



LEP200 since this machine will be operational in the relatively near future ||. We then 
turn to the proposed JLC/NLC/CLIC e + e~ colliders with possible center of mass energies of 
\/s = 500 GeV and 1 TeV |T0HT3|. We do not include any beamsstrahlung radiation effects 
in our calculation ||25|| . These effects are very much machine dependant (beam intensity, 
bunch geometry, etc . . . ) and known to be negligible at 200 GeV, and small at 500 GeV. 
However, although they can be quite important at 1000 GeV, there has been progress in 
strategies to minimize the effects of beamstrahlung radiation. 

II. CALCULATIONS AND RESULTS 

The diagrams contributing to the process e + e~ —>■ ^vivi are shown in Fig. 1. The main 
advantage of this process is that it depends only on the WW'y vertex. In addition, our signal 
(fig. 1(a)) should increase with energy, for two reasons: it is a t-channel process and should 
not decrease as fast with energy as the other contributions to the total process, especially 
when suitable kinematic cut are imposed to eliminate the on-shell Z contribution. Also, 
anomalous couplings, in general, become more important at higher energies. 

To evaluate the cross-sections and different distributions, we used the CALKUL helicity 
amplitude technique [f26[ to obtain expressions for the matrix elements and performed the 
phase space integration using Monte Carlo techniques [p7[| . The expressions for the helicity 
amplitudes are lengthy and unilluminating so we do not include them here. The interested 
reader can obtain them directly from the authors. To obtain numerical results we used the 
values a = 1/128, sin 2 9 = 0.23, M z = 91.187 GeV, Y z = 2.5 GeV, M w = 80.2 GeV, and 
r w = 2.1 GeV. 

The signal we are studying is an energetic 7 plus missing transverse momentum. The 
largest potential background is Bhabba scattering with a hard photon and the electron and 
positron going down the beam pipe. This should not be a serious problem given that we are 
interested in energetic photons and the luminosity monitors, based on Bhabba scattering, 
should be able to veto any such events. In order to take into account finite detector accep- 
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tance, we require that the photon be at least 10 degrees away from the beam line although in 
practice the cuts we use to enhance the signal are much more restrictive than this (typically 
~ 30°). 

In principle we should include QED radiative corrections from soft photon emission and 
the backgrounds due to a second photon either that is lost down the beam pipe or collinear 
to the hard photon being measured and therefore unresolved [23]. Although soft photon 



emission can reduce the cross section substantially, their inclusion does not substantially 
effect the bounds we obtain and therefore our conclusions. The effects of an unseen second 
photon turn out to be quite small. Since both these contributions depend on details of 
detector such as energy resolution and do not alter our conclusions we leave them out but 
stress that they must be included in detailed detector Monte Carlo simulations. 

The approach we followed was to examine various kinematic distributions, da/d cos d 1 , 
da/dEj, da/dEx-y, and kinematic cuts to find which ones optimized the sensitivity to anoma- 
lous couplings. In general, the tightest constraints were obtained by imposing cuts on the 
photon energy which eliminated contributions from on-shell Z° production (diagrams 1(d) 
and 1(e)). 

A. y/s = 200 GeV 

In Fig. 2 we show the distribution for yfs = 200 GeV for several values of k 7 and 
A T The most prominent feature of the distribution comes from the contribution of the 
on-shell Z° (diagrams 1(d) and 1(e)). From Fig. 2 it can be seen that the regions off the 
Z° resonance are most sensitive to anomalous couplings with the greatest sensitivity in the 
region above the Z° resonance. This is because that region probes the largest momentum 
transfer through the H^-boson t-channel propagators. We found that the tightest constraints 
could be obtained by imposing an angular cut of 35° < # 7 < 145° and on the two regions of 
£ 7 : 25 GeV< £ 7 <65 GeV and £ 7 >88 GeV. For 25 GeV< £ 7 <65 GeV the cross section 
is 0.14 pb which for an integrated luminosity of 500 pb -1 results in about a 12% statistical 
error. Similarly for E 1 >88 GeV a = 0.087 pb which gives a 50% statistical error. Monte 



Carlo studies of SLD type detectors give very crude estimates of systematic errors of 5% for 



cross section measurements | 29fl . Therefore the effects of including a 5% systematic error 
are not particularly important. The 68% C.L., 90% C.L., and 95% C.L. bounds that could 
be obtained with these cuts for y^i = 200 GeV and integrated luminosity of L=500 pb _1 
are shown in Fig. 3 with numerical values when varying one parameter at a time given in 
Table I. It is worth mentioning that the most efficient energy and angular cuts will vary 
slightly with the domain of k 7 and A 7 being probed. However, this dependance is mild and 
the limits will not vary much if one uses a fixed set of cuts. 

An important question for LEP200 is the effect on the physics reach of different center of 
mass energies. To gauge the change of the sensitivity to the TGV's for different energies we 
plot in Fig. 4 the 90% C.L. for yfs~ = 175 GeV, 200 GeV, and 230 GeV. Clearly, the higher 
the energy the tighter the constraints that can be obtained. Roughly speaking, increasing the 
cm. energy by 15% will increase the sensitivity to anomalous couplings by 30%. However, 
even for the highest possible energies at LEP200 of y/s = 230 GeV the obtainable bounds 
are not competitive with the iy-pair production process or for that matter with bounds 
obtained by the Tevatron experiments via associated Wj production, and will, at best, be 
a consistency check of other measurements. 

B. a/s = 500 GeV 

We next turn to an "NLC" type e + e~ collider with -^i = 500 GeV. We consider integrated 
luminosities of 10 and 50 /fa -1 . The photon energy distributions for a several values of k 7 
and A 7 are shown in Fig. 5. As before, the values of most sensitive to anomalous 
couplings are off the Z° resonance and k 7 and A 7 are sensitive to different values of E 1 . 
For example, k 7 is most sensitive to 50 < E 1 < 160 GeV while A 7 is most sensitive to 
177 < E 1 < 237 GeV. The bounds obtained from these observables are shown in Fig. 6 with 
numerical values obtained by varying one parameter at a time given in Table I. One could 
obtain additional information by using the transverse energy distribution of the photons 
with and without left-handed polarized electrons. For 50 < 2£ 7 < 160 GeV the cross section 

6 



is a — 0.27 pb which for an integrated luminosity of 50 fb 1 gives a statistical error of 
0.8%. Likewise for 177 < £ 7 < 237 GeV we obtain a = 0.024 pb and 5a stat /a = 3% and 
for E 1 > 245 GeV a = 0.0019 pb and 5a stat /a = 10%. Given the small statistical errors 
including a 5% systematic errors will have a considerable effect on the bounds that can 
be obtained. For example for 50 < E 7 < 160 GeV where the statistical error is smallest 
including a 5% systematic error reduces the sensitivity by roughly a factor of 5 while for 
177 < E 1 < 237 GeV where the statistical error is larger it has only a small effect and 
for E y > 245 GeV where the statistical error is largest the effect of the systematic error is 
neglegible. 

C.y/s = l TeV 

The final case we consider is a 1 TeV e + e~ collider. In Fig. 7 we show the photon energy 
distribution for several values of k 7 and A 7 . As in the 500 GeV case we find that k 7 and 
A 7 are sensitive to different values of E T For the photon energy range 60 < i? 7 < 220 GeV 
a = 0.39 pb which for an integrated luminosity of 200 fb~ x gives 5a stat /a = 0.3% while for 
245 < Ey < 490 GeV a = 0.05 pb 5a stat /a = 1%. Including a 5% systematic error weakens 
the bounds on the TGV's considerably, even more so at y/s = 1 TeV than at 500 GeV. 

In Fig. 8 we show the bounds that could be obtained using 60 < E 1 < 220 GeV and 
245 < E-y < 490 GeV and give some numerical values in Table I. The bounds obtainable on 
/t 7 are down to the percent level necessary to probe radiative corrections to the gauge boson 
vertices. Although the bounds on A 7 are also down to this level radiative corrections are 
about an order of magnitude smaller so that it is unlikely that deviations from the tree-level 
value could be observed unless there was a radical departure from the standard model by, 
for example, compositeness. 
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III. CONCLUSIONS 



We have examined the usefulness of the process e + e~ — > ^v\i>i for measuring the r )W + W~ 
vertex. The sensitivity of this process to anomalous coupling is greatly enhanced (by a 
factor of 5-8 generally) by eliminating the on-shell Z and by splitting the energy domain 
of the photon into more than one bin. The process turned out to be too insensitive at 
TRISTAN and LEP-100/SLC energies to obtain bounds competitive with recent results 
from the Tevatron. In fact, the same is also true for LEP-200 energies which at best will 
offer a consistency check of bounds extracted from H^-pair production at LEP-200 and the 
Tevatron. At higher energy e + e~ colliders, this process can lead to very stringent bounds, 
precise enough to test the TGV's at the level of radiative corrections. We used the high 
luminosities planned for at the high energy e + e~ colliders to estimate statistical errors. 
When we included reasonable estimates of systematic errors we found that the limiting 
factor in high precision measurements will likely be systematic errors not statistical errors. 
However, we based our analysis on optimizing cuts on and E^ and it is likely that 
performing a more detailed maximum likelihood fit of real data would make fuller use of the 
information in the data thereby obtaining tighter bounds. The challenge will be to reduce 
the systematic errors and one should be very careful with respect to the conclusions one 
makes by considering only statistical errors. 
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FIGURES 



FIG. 1. The Feynman diagrams contributing to the process e + e — > "fuv 

FIG. 2. The photon energy distribution, da/dE^ at y^s = 200 GeV. The solid line is for 
standard model values of k 7 and A 7 , the long dashed line is for k 7 = —3 and A 7 = 0; the dotted 
line is for k 7 = 5 and A 7 = 0; the dot-dashed line is for ac 7 = 1 and A 7 = —4; and the dot-dot-dashed 
line is for k 7 = 1.0 and A 7 = 4. 

FIG. 3. Sensitivities of the TGV's to anomalous couplings for y/s = 200 GeV and L=500 pb _1 . 
The horizontal curves are based on 25 < E 1 < 65 GeV and the vertically aligned oblongs are based 
on Et-j > 84 GeV. In both cases the solid curves represent 68% CL limits, the dashed curves 90% 
CL limits and the dot-dashed curves 95% CL limits. 

FIG. 4. Sensitivities of the TGV's to anomalous couplings at 90% CL for y/s = 230 GeV (solid 
curves), y/s = 200 GeV (dashed curves), and y/s = 175 GeV (dot-dashed curves) all for 500 pb _1 
integrated luminosity. The horizontal curves are based on 25 < E 1 < 80 GeV for y/s = 230 GeV, 
and 25 < E^ < 65 GeV for y/s = 200 GeV, while the vertical oblongs are based on E^ > 104 GeV 
for y/s~ = 230 GeV, E 1 > 88 GeV for y/s' = 200 GeV, and £ 7 > 72 GeV for y^ = 175 GeV. 

FIG. 5. The photon energy distribution, da/dE 1 for \/s = 500 GeV. The solid line is for 
standard model values of k 7 and A 7 , the long dashed line is for k 7 = 0.0 and A 7 = 0; the dotted 
line is for k 7 = 2.0 and A 7 = 0; the dot-dashed line is for k 7 = 1 and A 7 = —1.0; and the 
dot-dot-dashed line is for k 7 = 1.0 and A 7 = 1.0. 

FIG. 6. Sensitivities of the TGV's to anomalous couplings for y/s = 500 GeV and L=50 fb -1 . 
The horizontal curves are based on 50 < E 1 < 160 GeV; the vertical curves are based on 
E 1 > 245 GeV; and the upside-down U shaped curves are based on 177 < E 1 < 237 GeV. In 
all cases the solid lines are 68% CL. and the dot-dashed curves are 95% C.L.. 
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FIG. 7. The photon energy distribution, da/dEj for ^/s = 1000 GeV. The solid line is for 
standard model values of k 7 and A 7 , the long dashed line is for re 7 = 0.2 and A 7 = 0; the dotted 
line is for k 7 = 1.8 and A 7 = 0; the dot-dashed line is for k 7 = 1 and A 7 = —0.2; and the 
dot-dot-dashed line is for k 7 = 1.0 and A 7 = 0.2. 

FIG. 8. Sensitivities of the TGV's to anomalous couplings for yfs = 1 TeV and L=200 fb _1 . 
The horizontal curves are based on 60 < E 1 < 220 GeV; and the inverted U shaped curves are 
based on 245 < E 7 < 490 GeV. In both cases the solid lines are 68% C.L., the dashed lines are 
90% C.L., and the dot-dashed curves are 95% C.L.. 
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TABLES 

TABLE I. Sensitivities to k 7 and A 7 at 95% C.L. from the process e + e~ — > at a 200 GeV, 



500 GeV, and 1 TeV 


e+e colliders. 


The statistical error is based on the 


specified integrated 


luminosity. The entry 


denotes a bound too weak to be relevant. 




sfs = 200 GeV 


Observable 




5f a *(L=250 pb" 1 ) 


<5| tat (L=500 pb- 1 ) 


25 < £ 7 < 65 GeV 


<5k 7 


-2.5 


-1.9 


E 1 > 88 GeV 


<5k 7 


+2.3 
-3.2 


+1.9 
-2.6 




<5A 7 


+2.9 
-2.4 


+2.5 
-2.0 


E Tl > 84 GeV 




+2.9 
-3.3 


+2.5 
-2.8 






+2.8 
-2.4 


+2.4 
-1.9 


^ = 500 GeV 


Observable 




6( tat (L=10 ftr 1 ) 


<5|* a *(L=50 fb- 1 ) 


50 < £ 7 < 160 GeV 




+ .13 
-.12 


+ .055 
-.055 


177 < £ 7 < 237 GeV 


5k 7 


+ .16 


+ .08 
-.10 




<5A 7 


+ .16 
-.10 


+ .12 
-.06 


£ 7 > 245 GeV 


5X 7 


+ .17 
-.23 


+ .10 
-.17 


y/s = 1000 GeV 


Observable 




5f a *(L=50 fb- 1 ) 


<5| ia *(L=200 fb" 1 ) 


60 < £ 7 < 220 GeV 


5k^ 


+ .037 
-.036 


+ .018 
-.018 


245 < £ 7 < 490 GeV 


<5A 7 


+ .039 
-.015 


+ .034 
-.009 
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